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Anomalous Thermal Conductivity of Semi-Metallic Superconductors with 

Electron-Hole Compensation 

Hiroto ADACHI * and Manfred SIGRIST 

Theoretische Physik, ETH-Honggerberg, Zurich 8093, Switzerland 

The effects of low carrier density and carrier compensation on mixed-state thermal transport 
are investigated beyond the quasiclassical approximation. It is shown that, contrary to the usual 
observations, the interplay of the two effects leads to an increase in the thermal conductivity 
immediately below the upper critical field Hc2- Our result can account for the anomalous 
behavior of the mixed-state thermal conductivity near Hc2 recently observed in URu2Si2. 
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The discovery of the Nernst effect in the normal 
state region of underdoped cuprates^^ has initiated re- 
newed interest in magneto-transport phenomena in ex- 
otic superconductors. While a fluctuating vortex contri- 
bution may be a natural explanation for the unusual 
Nernst signal in cuprate, two other important aspects for 
magneto-transport phenomena were recently pointed out 
by Behnia and coworkers:^' the effects of electron-hole 
compensation^^ and low carrier density. '^^ The recent ex- 
periments on ultraclean URu2Si2 by Kasahara et aZ.^-* 
have motivated our study of these two points with regard 
to mixed-state thermal transport. The most intriguing 
result of ref. 4 is the increase in the low-temperature ther- 
mal conductivity below the upper critical field Hc2 , show- 
ing a hump structure.^) Usually in the low-temperature 
limit, the thermal conductivity of a clean type-II super- 
conductor decreases below Hc2j^ due to the enhance- 
ment of the Andreev scattering rate by vortices and the 
reduction of the density of states at the Fermi energy.^' ^' 

The key feature behind this anomalous behavior lies 
in the electronic structure of URu2Si2 introduced by the 
so-called hidden-order phase occurring below T* = 17.5 

10-12) jj-^ ^]-^jg phase, the carrier density is drastically 
reduced, as seen in several transport measurements. ^^"^^^ 
Furthermore, a nearly perfect i?^-dependence of the 
magneto-resistance without any sign of saturation sug- 
gests the compensation of electron and hole pockets of 
the Fermi surface. This interpretation is further sup- 
ported by the relatively small Hall angle. 

In this letter, we address the effects of the low car- 
rier density and carrier compensation on thermal trans- 
port in the mixed phase using a simple model of a two- 
dimensional s-wave superconductor. This allows us to ex- 
plain the unusual magnetic-field dependence of the ther- 
mal conductivity below Hc2 found in ref. 4 for URu2Si2, 
at least on a qualitative level. Since in this study it is 
necessary to go beyond the quasiclassical approxima- 
tion, we generalize the approximation scheme for obtain- 
ing the Brandt-Pesch-Tewordt Green's function^^^ valid 
near Hc2, by employing the formalism of Vavilov and Mi- 
neev^^'' originally developed to describe the mixed-state 



de Haas- van Alphen effect. 

We start by briefly reviewing the method of ref. 17. 
Gor'kov equations for a two-dimensional s-wave super- 
conductor under strong magnetic fields (k^ ~ c = h = 1) 
are as follows: 



^ Hoir) - u{r) - A(r)] G(r, r'; ie„) = S{r - r'), (1) 



where Ho + u + A^ (^"^+"^ J, G = ( ^t', g\ )' and 
En = 27rr(n+l/2) is the fermionic Matsubara frequency. 
The short-range impurity potential u{r) obeys the Gaus- 
sian ensemble u{r) = 0; u{r)u{r') = {l/m*T)S{r — r') 
with m* and r being the effective mass and the mean free 
time of quasiparticles, respectively. The single-particle 
Hamiltonian Tio is expressed as Hoir) = 2^Q'^ ~ A*' 
where Q = — iV + |e|A(r) with the vector potential 
A(r), and /i is the Fermi energy. For the magnetic field 
we assume the Landau gauge A{r) = Hxy, i.e., a uni- 
form field as justified near Hc2, if the Ginzburg-Landau 
parameter kql is large ( kql ^ 1 in URu2Si2^^''). For 
simplicity, we have dropped the Zeeman coupling term, 
although paramagnetic effects may be non-negligible in 
URu2Si2.^^^ The paramagnetic limiting effects are, how- 
ever, beyond the scope of this letter. 

The single-particle Hamiltonian Hoir) can be diago- 
nalized with the eigenvalues Cat = uJc{N + 1/2) — /i and 
the eigenfunctions^°^ 
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where LOr 



\e\H/m*, (Pn{x) ^ 



X + ( 



-(y-A^fe) + igay 



A 



(2) 



y/2'^N\^ 

A = (|e|i/)~^/^, Hjm{x) is the N*'* Hermite polynomial, 
and q represents the quasi-momentum in the magnetic 
sublattices.^^^ We have chosen a rectangular latticc"'^^-' 
with edges ax = a and ay = 2a {a — \/ttX) ■ We introduce 
the magnetic sublattice representation 

G'(ri,r2;ie„) 
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= yi (pNi{ri\qi)GNuNAQi,Q2;i£n)(l)*N^{.r2\q2), 
— I <ii 

"1=0 12 
"2=0 

(3) 

with the shorthand notation f = f^^"', ^ f'^^^/o 

Jq J —IX j a liT J— 7r/2a 27r 

As in ref. 17, wc assume a square vortex lattice, and 
neglect, for the moment, the impurity potential for sim- 
plicity. 

Using the magnetic sublattice representation, the nor- 
mal component of G can be obtained from eq. (1) as 
follows: 

<^JVi,JV2(gi,Q2;ien) = G%^{\£n){'iT^f5{qi,q2)5NuN2 

oo 

-G%^{ien) 

^N^,N4iqi;i£n) 

N4=0 



of impurities on C?^(i£„) is included by the replace- 
ment ie„ ^ lEn - J2ns=o GN.M^n), if wc ucglcct the 
anomalous self-energy^^^ assuming the ultraclean limit. 

Inserting the approximate single-particle Green's func- 
tion into the Kubo formula, the thermal conductivity 
of a clean type-II superconductor in the low temperature 
limit is given by^' 



T 



dCNGf^^^'^G 



jv-i-1 Y 



fp/A 



-) 



)]}. (8) 



X GN^^N2{qijq2',i£n), 



(4) 



uf/A ^ i^f/A 

A/\/2, t;F is the Fermi velocity, and I{z) = 
J^^dt^—^. The retarded (R) and advanced (^4) 
Green's functions at zero energy are given by 



where A 
1 



'^Ni,N4iqi;^£n) = ^ '^Ni,N3{qi)^*N4,N3{qi)G%^{-i£n 

N3=0 

(5) 

where G%{iSn) = (i^n— Ca^)"^ is the normal state Green's 
function, and A^^^Nziq) given by^'') 
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where G^^(ie„) = (GO,^(i£„)-i - 
S^^(i£„) = ^5 



-|A|^ 



(jVi-A^3)' 
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Based on ref. 16, wc replace Sat^^jv^ (qi; i£„) in eq. (4) 
by its average over qi, because AjVi,Ar3(q'i) appears here 
in a gauge-invariant manner as |Ajvi,jv3(qi)P and the 
<7i-dcpendcncc has a small effect near Hc2- Note that 
the qi-average in the magnetic sublattice representation 
corresponds to the spatial average in the real-space rep- 
resentation. Now, we obtain 

GjVi,jV2(Qi,92;i£n) = {2Trf6{qi - q2)^JVi,JV2G^i (i^n), 

(7) 

(ie«)) , 

and we have used the Gaussian approximation^^^ in 
eq. (6) because we are interested in the situation A^f = 

^J,/u!c > 1. 

The Green's function (7) is a generalization of the 
Brandt-Pesch-Tewordt Green's function^^^ to a case be- 
yond the quasiclassical approximation, and it allows us 
to study the effect of the non-zero cyclotron frequency 
iVc- In our approximation, F and A include the off- 
diagonal elements in the Landau level space spanned by 
eq. (2). This is important, since the off-diagonal matrix 
elements incorporate the phase coherence of the Cooper 
pairs. In the limit of small lOcT, our approximation re- 
covers the Green's function of ref. 16 and the iJc2-line 
of ref. 23. This is in contrast to other approximations 
that lead to a reentrant behavior for i?c2-^^'^^^ In this 
respect, wc believe that our treatment is adequate for 
the qualitative description of URu2Si2. Finally, the effect 



The sum over Landau level index N is replaced by the 
integral over ^jv, which is obtained by transforming the 
A'^-summation by the Poisson summation formula, and 
only the non-oscillatory terms are retained. The density 
of states at the Fermi surface, g^v, is determined self- 
consistently."'^^'' 

We now consider a compensated superconductor by in- 

troducing the pairing interaction Vpair = ( yc- -)', y(h,h) )> 

where V^'^'^^ denotes the pair scattering amplitude be- 
tween bands a. and /?. Neglecting the intcrband impu- 
rity scattering, the thermal conductivity of such a su- 



-(e) 



where the 

''J 



perconductor is given by Ki 

contributions of the electrons k\^^ is given by eq. (8) 
and that of the holes k^^^ is given by eq. (8) after re- 



placing c 



N 



-(n in the integrand. We fix the length 



(e) 



scale to the coherence length of the electron band, 

/2'kTc^\ assuming that the electron band dominates 

siiperconductivity.^-* Then, wc specify the following pa- 
rameters: kp^^Q^^ is proportional to the area of the elec- 
tron Fermi surface; fcp'^V^F ' = ^/nfl^^JnJ^ measures the 
degree of compensation; /i^p/Co**^ and ^imp/^imp introduce 
the mean free paths of both carriers; Ag'^^Ag''^ (or equiv- 
alently Vpair) denotes the ratio of zero-field gaps. We keep 
■j/(e,h)^ y(h,e) ^ which rcsults in the field dependence 

A(-/^){H) = A^^'''^^l-H/H,2 near H,2. 

In the comparison between the theoretical and exper- 
imental results for thermal conductivity in a supercon- 
ductor with moderate strength of WcT, it is important 
to notice that experimentally the thermal resistivity ten- 
sor is measured, and not the thermal conductivity 
tensor k?"^^ Thus, the experimental longitudinal thermal 
conductivity ^ is given by 



- K -I- 



(10) 



where we used the relations Kyy = Kxx and Kyx = —Kxy 
First, we address the case of perfect compensation in 
the normal and superconducting states, i.e., the states 
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Fig. 1 . Magnetic- field dependences of thermal conductivity. Here, 



fe+h) 



(h),(h) 



) is the normal state thermal 



conductivity at zero magnetic field. V^air = V^'^'"^ 
used. 



Fig. 3. Main panel; Calculated thermal conductivity for esti- 
mated parameters of URu2Si2. T4>air = V^"''^') (o 35 "o*)" 
set: Magneto- resistance as a function of H"^ for several values of 
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Fig. 2. Enhancement of the thermal conductivity by the devi- 
ation from perfect compensation in the superconducting state. 



is depicted as a function of magnetic field for different ra- 
tios of the two gaps. With a decreasing ratio Aq^^/Aq''\ 
the increase in k^J^ below Hp2 is enhanced. Evidently, the 
violation of compensation in the superconducting phase 
enhances the hump structure below Hc2- 

Finally, we compare our results with the mixed-state 
thermal conductivity of URu2Si2. Note that our anal- 
ysis yields only a qualitative understanding because 
we employ a simple quasiparticlc picture with a two- 
dimensional isotropic Fermi surface and neglect the para- 
magnetic effect as well as the correlation effects. We fix 
the parameters as follows. From the values Hc2 ~ 2.8 
T for if II c together with A;f ~ 1-1 nm"^,*^) and taking 
into account the fact that there is a substantial paramag- 
netic effect -"^^^ in this material, wc have a rough estimate 
fcp^o — 6.5. For the mean free path, we use ^[^p/Co'^'' = 60 



top to bottom. 



in which the relations /k^'^^ 



1 l^^^ 

^1 '■imD/''ii 



1, and 



' imp/ imp 

1 hold. Figure 1 shows the magnetic-field 

(c)e(c) 



dependences of for several values of k^^^^' S,)^' ■ For a 
large value of k'^^'S^^^^ (= 10), the calculated de- 
creases for fields below Hc2- This is consistent with the 
well-known behavior^^^^ found in quasiclassical calcula- 
tions {k^^S,^^ oo). However, as the carrier density (or 
k^^Q^^) is reduced, a new property appears. On lower- 
ing the magnetic field, initially increasc;s below Hc2 
and then decreases, forming a hump structure. In case 
of a single carrier-type without electron-hole compensa- 
tion, does not show this kind of hump structure even 
for small values of k^^l^^ because of the second term in 
eq. (10). 

The hump structure of below Hc2 becomes even 
more pronounced if the sizes of the energy gaps are dif- 
ferent, i.e., A[,''V^o°^ 1- In this case, the second term 
in eq. (10) becomes nonzero below Hc2 and enhances the 
size of the hump. This is actually seen in Fig. 2 where nf^ 



imp/ -imp in order to reproduce the observed 

magneto-resistance (Ap3,a,(10T)/p^^(0T) w 300) with a 
nearly perfect H^-dependence (see the inset of Fig. 3). 
This choice of parameters is consistent with the sce- 
nario^-* in which the Hall effect is dominated by the 
light hole band. Further, the positive Hall coefficient 
and the measured value-^^^ pxy{lOT^) / pxx{OT^) ~ 50 gives 
an estimate fcp^V^F*^ = 1.01. Finally, we have assumed 
that a larger gap is formed in the electron band follow- 
ing the discussion in ref. 4. For the gap ratio, we use 



A (h) / A (e) 
^0 /^O 



0.31. 



The main panel of Fig. 3 shows the calculated thermal 
conductivity of URu2Si2 as a function of the magnetic 
field. The peculiar magnetic-field dependence of the ther- 
mal conductivity with the hump structure below Hc2 is 
reproduced, with its size being slightly larger than the 
measured one in our calculation.^' The neglected inter- 
band impurity scattering would reduce the size of the 
hump. It should be noted that while the size of the hump 
is modified by changing the parameters, the appearance 
of the hump itself is robust for parameters reproducing 
the magneto-resistance data.^' In the experiment, the 
thermal conductivity increases sharply below Hc2, which 
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(a) Impurity scattering 



(b) Impurity scattering 




drift motion 

Andreev scattering 
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Fig. 4. Schematics of the quasiparticle trajectory where Rc = 
A^fcp is the radius of the cyclotron motion of quasiparticles; (a) 
high field {wcT > 1), (b) low field {wcT < 1). 



most probably results from a paramagnetic limiting ef- 
fect in this material, which we neglect here. In our model 
case, the hump develops more gently. To treat these other 
effects is beyond the scope of our study and therefore we 
leave it for future studies. 

The special scattering features provide some physi- 
cal insight into this behavior. For the superclean limit 
LOcT ^ l,^^-* the heat transport along — VT by quasipar- 
ticles is only possible through scattering (see Fig. 4(a)), 
as otherwise the quasiparticles would experience only the 
drift motion perpendicular to — VT. In the normal state 
at low temperatures, this corresponds mainly to impu- 
rity scattering. In the superconducting mixed phase, an 
additional contribution is derived from Andreev scat- 
tering.*' Hence, at high magnetic fields, the Andreev 
scattering tends to increase the thermal conductivity be- 
low Hc2- At lower fields (smaller lOct), the Andreev scat- 
tering plays a different role. The particle and hole tra- 
jectories tend to retrace (Fig. 4(b)). Because both quasi- 
particle types carry heat, a compensation occurs, which 
decreases the heat current parallel to — VT and nf^ de- 
creases at lower fields. The overall behavior leads to the 
characteristic hump feature. Note that this picture of the 
single-carrier Kxx is applicable to k^-^ because electron- 
hole compensation leads to the suppression of the second 
term in eq. (10). In the ordinary case without the com- 
pensation, this picture cannot be applied to nf^ because 
the second term in eq. (10) makes an important contri- 
bution to Avf-^ at high fields, veiling the hump below Hc2- 

In conclusion, we have examined the effects of the low 
carrier density and the electron-hole compensation on 
mixed-state thermal transport. The interplay of these 
two effects leads to the appearance of a hump structure of 
Kxx{H) below Hc2, as observed in URu2Si2, whose elec- 
tronic states incorporate the two features in the hidden- 
order phase. Our study provides a natural explanation 
for the unusual magnetic- field dependence of Hxx{H) 
near Hc2- Moreover, it demonstrates that URu2Si2 pro- 
vides a rare chance to examine the physics of the mixed 
phase in the superclean limit, resulting in intriguing 
magneto-transport phenomena. 
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